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I. INTRODUCTION 



The form factors of heavy-to-light transitions at large and intermediate energies are 
among the most important applications of QCD light-cone sum rule (LCSR), since the 
validity of the LCSR approach is restricted to the large meson energy {Ep >> Aqcd) via 
the relation = m\ — 2mBEp. In literature there are several approaches to calculate the 
B —>■ light meson transition form factors in addition to the QCD LCSR approach, such as 
the lattice QCD technique and the perturbative QCD (PQCD) approach. These approaches 
are complementary to each other, since they are adaptable in different energy regions, and 
by combining the results from these three methods, one may obtain a full understanding of 
the B — light meson transition form factors in its whole physical region 

mm 

. Since 

the LCSR is restricted to small and moderate g^, a better LCSR shall present a better con- 
nection to both the PQCD and the lattice QCD results, and then a better understanding of 
these form factors. 

How to "design" a proper correlator for these heavy-to-light form factors is a tricky 
problem. If the correlator is chosen properly, one can simplify the LCSR greatly. As for the 
B — > light pseudoscalar mesons, the commonly adopted correlators are usually defined as 

nj(p,g) = t J d^xe'''■^{P{p)\T{q{xh,b{x),m^mhl5q'm\0) (1) 



and 



Kip,q)=i J rfW-"(P(p)|r{g(x)«cr^,g'^6(a;),K0)^m,75g'(0)}|0), (2) 



where q{x) and q'{0) stand for the light quark fields that form the pseudo-scalar mesons. 
By taking such conventional correlation functions, it has been found that the main uncer- 
tainties in estimation of the B P form factors come from the different twist structures 
of pion/kaon wave functions, and most irnportantly, the twist-2 and twist-3 contributions 

nriQ 

should be treated on the equal footing p, la, 17|]. Thus one has to calculate both the twist-2 
and twist-3 contributions up to one-loop accuracy in order to obtain a consistent one-loop 
estimation of the form factors. 

On the other hand, by taking proper chiral currents into the correlator, one can directly 
eliminate the most uncertain twist-3 terms, and then only needs to calculate the twist-2 
o„„.ibu«„„ to „ext.to..ead„, o.de. (NLO) acc.acy BflQ. M p..,., *e vecto. 
form factors fB^T^Kil"^) has been calculated with the chiral current in the LCSR up to 



NLO 



9|. It can be found that the scalar and penguin form factors fs'^nKil'^) shall be 



important in due cases, e.g. the penguin form factors shall give sizable contributions to 



B Pl^l^ or B —>■ K*'y So it is interesting to extend the previous study to all the 
B ^ P {P stands tt, K and all light pseudoscalar mesons) transition form factors f^'^'pi.'f') 
with the chiral current in the LCSR up to one-loop accuracy. Furthermore, it maybe also 
interesting to know to what degree the different choices of correlator shall affect the final 
LCSRs, which is another purpose of present paper. 

The paper is organized as follows. In Sec. II, we present the calculation technology to 
obtain the LCSRs for the B ^ P transition form factors /^ijp(g^) with chiral currents, 
where the S'f//(3)-breaking effects for the kaonic case will be explained in due places. Nu- 
merical results and discussions are presented in Sec. Ill, where the uncertainties of form 
factors under the present LCSRs shall be discussed. The comparison with other approaches 
will be presented in Sec. IV. Sec.V is reserved for a summary. 

II. CALCULATION TECHNOLOGY FOR THE B ^ P TRANSITION FORM 
FACTORS WITH PROPER CHIRAL CURRENTS 

A. A definition of /^^'p 

Based on the previous calculation about the transition form factor B —>■ n/K, we present 
the formulae for the B —>■ P transition form factors for generality such that these formulae 
can also be conveniently extended for other light pseudo-scalar form factors like B ^ rj and 
B ^ 1]' form factors. With default, we adopt the chiral limit pp = mp = 0, but point out 
how to include the S'f//(3)-breaking effects for the B ^ K form factors in due places, i.e. 
the dominant 5'[//(3)-breaking effects will be discussed with the newly obtained K meson 
distribution amplitudes The hadronic matrix elements for the B ^ P transition form 
factors are parameterized as 

= 2f^^p{qVp + [f^-.p{q') + fE^p{q')]q, , 

{P{pp)W^a,,q''b\BipB)) = /L^i^l [p.qq^- q^pj (3) 

mp -|- mp L -I 
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with P representing the pseudoscalar, = (ps + Pp)^, q^i = {pb — pp)fi, and fB^pi.'f')) 
fB~*p{Q^)i fB^pio^) stand for the vector, scalar and tensor form factors respectively. It can 
be found that the scalar form factor fB^pi'f) satisfies the following relation: 



fi^p^e) = f^^pig') + ^^T^fB^pi^')- (4) 

lll^Q 'lip 

As for the LCSR calculation, different to the conventional choice of the correlation func- 
tions as shown in Eqs. f|T][2l) . we choose the following chiral currents in the correlation func- 
tions, 
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nj(p,g) =ij rf^xe*^-^(P(p)|T{g(x)7M(l + 75)6(^)J(0)^m,(l + 75)g'(0)}|0), 

= n+(g2, (p + qf)p^ + U-{q\ {p + qf)q^, (5) 



= If{q\ ip + qf) [(P ■ q)q, - q^P,\ , (6) 
where P = p + 2q. 

We calculate the form factors f^'^pi^f) following the same calculation technology as 



described in Refs. 



9|, where the vector form factors f^^^K^'f) have been calculated. For 
such purpose, we first give a simple extension to fs-^piq"^) in the large space- like momentum 
regions {p + q)"^ — rnl <^ and ^ ml for the momentum transfer, which correspond to 
the small light-cone distance ~ and are required by the validity of OPE. And then, we 
present the newly obtained results for the scalar and tensor form factors. 

B. A simple extension to f^^p within LCSR 

The vacuum-to-meson matrix elements in terms of the pseudo-scalar's LC DAs of different 
twist can be expanded by contracting the 6-quark fields with the help of the full 6-quark 
propagator within the background field: 

(0|T6(x)6(0)|0) = ^ f -^^e-^'^^^, - igs I T^^e"^^^ • 



1 

dv 







(27r)4 k^-ml (27r)4- 

^ (7) 



where only the free propagator and the one-gluon terms are retained, G^^ stands for the 
background gluonic field strength, and Qs denotes the strong coupling constant. The in- 
variant amplitudes n"*" can be obtained by substituting the 6-quark propagator and the 
corresponding LC wave functions, and completing the integrations over x and k. 



The OPE results for the invariant amphtudes n"*" can be represented as a sum of LO and 
NLO parts: 



An 



(8) 



where nQ'(g^, (p + g)^) and 'n.i{q^, {p + g)^) stands for the LO and the NLO contributions 
respectively. As for the LO invariant amplitude, we obtain: 



U+{q',{p + qf) = 2fpml 



1 du (fp{u) 



u 



A 



1 du 



mt 



u 



3 2A3 



1 du 



;Gip(u) 



+ 



1 f ^ 2^4p(ai) + 2^4p(a 
dv / Dai 



mA^ 
$4p(ai) - ^4p(ai) 



A2(ai + f tts)^ 



(9) 



where the parameters are defined as: A = s — {p + g)^ (s = [g^ + (m^ — g^)/-u]), G'4p(m) = 
— Jq dvtp4^p{v) and Dctj = (iai(ia2(ia35(l — ai — 0:2 — 0:3). Here ifp is the twist-2 LC wave 
function, and 4>p{u), ip4p{u), \l/4p(aj), \^4p(aj), $4p(aj) and <l>4p(aj) are twist-4 LC wave 
functions defined in a same way as the pionic case that have been defined in Ref.fl2|, whose 
explicit forms are put in the APPENDIX A. It is found that only the twist-2 and twist- 
4 contributions are contained in the above expressions, and the twist-3 terms are rightly 
eliminated by taking the present adopted chiral currents within the correlators. 

Since the most uncertain twist-3 contributions are eliminated and the twist-4 contribution 
itself is quite small, so we only need to consider the NLO correction to the twist-2 terms. 
The NLO invariant amplitude n^*" for the twist-2 contribution can be written in the following 
factorized form: 



{P + qf) = -fp t duTt{q\ {p + g)^ u)Mu) 
Jo 



(10) 



where by taking nib to be the b-quark pole mass, the NLO hard scattering amplitudes Tf*" 
can be written as 
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u 



+ 



1-p 
4 



^2 -In 4^ 



P 



[2G{p) - Gin) - G(r2)] 



(ri - r2 



^ [Gip) - G(rO] + ^—^ [G{p) - G{r,)] 



u 



u 



+ 2 



p+ (1 -p)ln(l -p) 4 (l-r2)ln(l-r2) 



p-r2 



p^ 1 - P ^2 

;i-p)ln(l-p) (l-r2)ln(l-r2) 



P 



r2 



(11) 
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with 



u 



rx 

1 - u, p = [ri + u(r2 - ri) - u(l - u)Mp/ml], Li2(x) = - — ln(l - t), 

Jo t 



mt 



G (p) = hh{p) + ln^(l - p) - ln(l - p) ( 1 - In -| 1 , 



(12) 



where the dilogarithm function Li2(x) — — Jq y ^^(1 ~ '^i = ?^/"^& ^-iid ^2 = (p + qT/ttT'I- 
Next, the QCD LCSR for fB^p{q^) can be schematically written as 

1 



«0 



(m|-s)/M2 



gV'"s 



(13) 



where Pt2(^)?^) ^l^® contribution from the twist-2 DA and PtaIq^) is for twist-4 DA, Jb 
is the B-meson decay constant. The Borel parameter and the continuum threshold Sq 
are determined such that the resulting form factor does not depend too much on the precise 
values of these parameters; in addition the continuum contribution, which is the part of the 

dispersive integral from sq to 00 that has been subtracted from both sides of the equation, 
should not be too large, e.g. less than 30% of the total dispersive integral. 
As for the LO twist-2 and twist-4 contributions, we obtain 



fBf^^p{q')\LO 



mil 



ml 



p "'s 

e 



due~ 



u 



Jo J {a^ + va^YM^ 



2*4p(ai) + 2^ip{ai) - $4p(ai) - $4p(ai 



(14) 



where A = for Mp = 0; A = V(^o-^--M|)-+4M^».g-.-)-(.o-.--M^) 

As for the NLO twist-2 contribution, it is convenient to write the NLO Pt2{^iQ^) the 
following form 

fp [asCF\ /"^ 



where 



Pt2\s,(1 )\nlo 



TT V 47r 



/ du(pp(u, p)lmT'^ 
Jo 



1 1 —2' ^2'^' A'' ) 



mt mt 



— Im,r+ = 9{1 - p) 
+ 0{p-l) 



^2(^2) 



+ 



1 — ri 



p-1 

L2(r2)-2Li(p) 



{r-2 - ri)(r2 - p) 



P 



+ 



1-ri 



(r2 - ri)(r2 - p) 
1 2 ^ 



r2 - 1 
(r2 - p)r2 

Li{r2) 



ri-p r2-p 2p V p r2 



(15) 



In 



r2 - 1 
1 — ri 



r2 



ri 



m 



71 



- 2^^(^7^j +Li2(ri)-3Li2(l-r2) + 3-y 
for the case of ri < 1 and r2 > 1. The operation " + " is defined by 



dpfip) 



1-p 



dp[f{p) - /(I)] 



1-p 



The two functions Li[x) = In 



2 

b_ 

X /i^ 



(x-i) 



(x-i) 



2 

o 



(16) 



(17) 



- are introduced 



1 and L2{x) = In 

to make the formulae short. The above formulae are derived in the Feynman gauge and 
by regularizing both the ultraviolet and coUinear divergences by the standard dimensional 



regularization in the MS scheme. 



C. Calculation of /°'^p within LCSR 

For convenience, we calculate the combined function fs^pi.'f') = fs^pi.'f') + fB~>pil 
first and then derive f%^p with the help of Eq.(jl]). The OPE results for the needed invariant 
amplitudes n~'^ can be represented as a sum of LO and NLO parts: 

^-'^(g^ (p + qf) = Ho '^(g^ {p + qf) + ^^^^(g^ (p + q?). (18) 

where IIo '"^(g^, {p + qY) and n]~'"^(g^, {p + qY) stand for the LO and NLO contributions 
respectively. As for the LO invariant amplitudes, we obtain: 

^o(g^(p + g)^) = 2fpml / ^-r-,G,p{u), (19) 

Jo 

Hl{q\{p + qf) = 2m,fp 



^du(pp(u) 1 / 2m? \ , , , 



u A Jo Au^A? \ uA ^ 
+ t dv /•2)«.2^4p(ai)-(l-2t;)$4p(a^) + 2(l-2t;)^4p(ai)-$4p(ai 



/o J A2(ai + m3)2 J^^°^ 

Similar to the case of no'(g^, {p + qY)-, one may also observe that only the twist-2 and 
twist-4 contributions are contained in the above expressions, and the twist-3 terms are 
rightly eliminated by taking the present adopted chiral currents within the correlators. The 
NLO invariant amplitude LIj^ ' for the twist-2 contribution can be written in the following 
factorized form: 

^^^(g^(p + g)') = -fp C duT^'^{q\{p + q)\u)vp{u), (21) 

Jo 



where by taking mt, to be the b-quark pole mass, we have 

2(ri - r2)[n + (1 - ri) ln(l - n)] 2(1 - n){n + r2) ln(l - n) 



Ti (ri,r2,n) 



r2(l - p) 



rKn - p) 



^ 4(l-r2)ln(l-r2) _ 2(1 - p)(r2 + p) ln(l - p) ^ 2{n - r^) ^^2) 



r2{p - r2) 



u{p - r2)p2 



rip 



and 



mt 



1-p 
4 



3-21n-| 



1-P 



[2G(p) - G{n) - G{r2)] 



(ri - ra)^ 



^ ''[G(ri)-G(p)] + ^[G(r2)-G(p)] 



u 



u 



^ ^ ln(l — r2) ^ ln(l — ri) 



1 - p V r2 



ri 



+ 4 



'-^^ \n ("i^ V 4 ("-1^ Vn ("i^" 



7i(p-n)/ \l-p) \ip-r2)r2j \l-p 
- 2 ( P±ML^] + f -4n + 2r2n+4r2\ 
\ P J \ nr2p J 

Schematically, the QCD LCSRs for f^^p can be written as 



1 r^o 



fsfB^piq') = -2 L e^™^ 



rriB + mp p ^(^2^_,)/a./2 
2m% 



pT2{s,q ) +Pt4(9 ) , 



pT2is,q )+ pTiiq ) 



(is 



(23) 

(24) 
(25) 



As for the LO twist-2 and twist-4 contributions, with the help of the Eqs. fll9ll20p . we 
obtain 

G4p(u) 



fBfB^p{q')\LO 
fBfE^p{q')\LO 



2-/2 2 ^ 

—u(q —um,p) 

e'M^ I due i^A?5 



2ml fp ™^ '"^ 



m 



B 



A 



(26) 



{ms + mp)mjj,^<ij^^^^^_ 



ml 



2-/2 2 \ 

1^ —u{q — um p) 



Vp{u) _ 04p(w) L ml \ 



u 



+ 



J («! + fasj^M^ 



X 



2^4p(«i) - (1 - 2^;)$4.(aO + 2(1 - 2t;)^4p(«0 - '^'4p(ai)J |. (27) 
From the above equations, we immediately obtain the relations among 

ftipiq'^) at the LO 

and up to the twist-3 accuracy, i.e. 



f-B^p{q') = -n^p{q') and /J_.p(g^) = "^^^^ft^Aq") 



mb 



(28) 
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which agree with the conclusions drawn in Ref. 13|. Moreover, with the help of Eqs. (Hll28l) . 
we obtain 



I B^P 



m 



B 



m, 



As for the NLO twist-2 contribution, the NLO Pj^^s^q^) can be written as 



Pt2 {s,q )\nlo = -— (^^^ 



where 



— InisTi 



+ 



0(1 -P) 

2(r2 - 1 
r2 



and 



Air 



lm,T^ = 9{1 - p) 



+ 0{p-l] 



2(1 -ra) 
+ 5{p-l) 



L2{r2) 



ri - p 

^2 



1 - 



ri - r2 _ (2 - r2)(r2 - ri) 
(ri - l)(ri - ln(l - ri) 



ri 



1 — ri 



p-1 
L2(r2) - 2Li(p) 



(r2 - ri)(p - r2) 



i:i(r2) 



^^2 



p-1 
^I(P) 



1 — ri 



'^2(p - 



ri + r2 - p - I 
(ri -p)(r2 -p) 

r2 - 1 



+ (r2 - ri)(p - r2) 
ri - 1 2(r2 - ri) + rir2 



+ 5(p-l 



- 1- 



In 



ri 



ri(p - ri) 
(^2 - 1)^ 



2rir2p 



+ 



2p2 



1 \ TT^ 

In^^^ ^ (lnr2 + 1 +3 - — 

1 - ri V r2 / 2 



ln(l - ri) - 2 In 



+ Li2(ri) - 3Li2(l - 



(29) 



/ du(l)p(u,p)lmT^''^ ( ■^,^,u,p] , (30) 
JO mi ) 



(31) 



(32) 



ri / \ p^ 

for the case of ri < 1 and r2 > 1. 

As a cross check of the above NLO formulae for the twist-2 contributions, it can be found 



that our present results for /^I^^'J agree with Ref.|l4] by transforming the formulae for the 



MS 6-quark mass to be the ones for the 6-quark one-loop pole mass, except for an overall 
factor 2 ^ . 

Here similar to the treatment of Refs. N, QJ^, M , we have adopted the 6-quark pole mass to 
do the calculation. Refs. [M, [ITJ have argued to use the 6-quark MS running mass other than 
the pole mass. Numerically, we shall show in due places that if properly choosing the possible 
ranges for the undetermined parameters, these two treatments are in fact equivalent to each 



The overall factor 2 comes from the different choices of correlation function. 
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_ 


so (GeV2) 


M2 


fB (GeV) 


ruh = 4.75 (GeV) 


33.0 


2.48 


0.192 


mb = 4.80 (GeV) 


32.6 


2.28 


0.169 


ruh = 4.85 (GeV) 


32.3 


2.10 


0.146 



TABLE I: The value of Jb (in units: GeV) within the LCSRs with chiral currents up to NLO, the 
corresponding formulae can be found in Ref.^], where mi, is taken to be the 6-quark pole mass. 

other within reasonable uncertainties. We prefer to take the pole quark mass, sine the pole 
quark mass is universal that can be determined through proper potential model analysis 
or through lattice QCD calculation, while the running quark mass is process dependent, 
i.e. depends on the renormalization scheme and the renormalization scale of a particular 
process. 



III. NUMERICAL RESULTS FOR fefi'J^Ki'l^) WITHIN THE QCD LCSR WITH 
CHIRAL CURRENT 

A. Parameters and distribution amplitudes of the light mesons 

First, we specify the input parameters used in the LCSRs for B ^ it and B ^ K 
transition form factors. For the needed meson masses and the light mesons' decay constants, 
we adopt the center values as listed by the Particle Data Group [18 1 



U = 130.4MeV, fK = 155.5MeV, 

Mb = 5.279GeV, = 139.570MeV, Mk = 493.667MeV. 

As has been argued in the last section, we shall adopt the 6-quark pole mass to do 
numerical calculation throughout the paper. As for the value of fs, to be consistent with 
the present calculation technology, they should be determined by using the two-poiiit sum 



rule with proper chiral currents up to NLO. Such a calculation has been done in Ref. [s], the 
interesting reader may turn to Ref.|3i for more calculation detail, and here we only quote 
some typical results as shown in TabUl where the one-loop pole mass rrib is taken to be 



(4.80 ±0.05) GeV (l9|. 



Naively, the leading twist-2 DAs 0^ and (px can be expanded as Gegenbauer polynomials 

10 



0.3 




a;(1GeV) 



FIG. 1: 02 (IGeV) and aJ(lGeV) as determined from the two constraints adopted in the body of 
the text, where the rhomboid stands for the allowable range. 

as shown in APPENDIX A. The first two Gegenbauer moments, e.g. 03 and for pion and 
af and af for kaon, have been studied with various processes. We adopt two constraints 
for a^(lGeV) and aJ(lGeVL e^. a^(lGeV) +<(lGeV) = 0.1 ±0.1 ^ and -|a^(lGeV) + 



T|aJ(lGeV) + | = 1.2±0.3 [5|,|2l|], such that the allowed values of and are correlated and 



given by the rhomboid shown in Fig. ([T]) . Note here we do not adopt the wider range of 02 = 
0.25±0.15 as suggested by Ref.[l2|, since we prefer a more asymptotic-like pion DA as favored 



by a very recent QCD LCSR analysis of B ^ n vector form factor [27. 



moment af has been studied by several references, e.g. Refs. 12|, 22I, 3, 24, 25, 26| and etc. 
For convenience, we quote the values for the twist-2 Gegenbauer moments of kaon as obtained 
from the average of those obtained in literature to do the discussion, af- (IGeV) = 0.06±0.03 
and af (IGeV) = 0.25 ± 0.15 (l2j. 

Furthermore, for the twist-2 DAs, we do not adopt the Gegenbauer expansion flAip . since 
its higher Gegenbauer moments are still determined with large errors whose contributions 
may not be too small, i.e. their contributions are comparable to that of higher twist struc- 
tures |3i]. As a compensation, we adopt the suggestion of deriving the pion and kaon DAs 
from their corresponding WFs by integrating over the transverse momentum j^. And the 
twist-2 pion and kaon WFs can be constructed on their first two Gegenbauer moments and 
on the BHL prescription 28|], i.e. 



le first Gegenbauer 



^^{x, kx) = [1 + B^C^^^\2x - 1) + C^Ct'\2x - 1)] 



3/2/ 



x(l — x) 



exp 



' k5_ + 
x{X — x) 



2\ n 



(33) 
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«2(Mo) 


0.00 


0.115 


0.230 


04(^0) 


0.00 


0.092 


-0.015 


-0.120 


-0.030 


A^{GeV-^) 


226.0 


196.7 


199.4 


199.1 


173.8 




-0.079 


-0.024 


-0.018 


-0.014 


0.043 




0.027 


0.073 


0.012 


-0.050 


-0.00656 




0.902 


0.862 


0.868 


0.870 


0.832 



TABLE II: Pion twist-2 wavefunction parameters for some typical Gegenbauer moments, where 
HQ = IGeV. Note the obtained WF parameters are for /i = 2.2GeV. 





0.09 


0.06 


0.03 


af (/io) 


0.40 


0.25 


0.10 


0.40 


0.25 


0.10 


0.40 


0.25 


0.10 


AKiGeV-^) 


171.2 


209.3 


253.8 


172.6 


211.8 


255.9 


173.9 


213.5 


258.1 


Bk 


0.0845 


0.0732 


0.0588 


0.107 


0.0966 


0.0825 


0.130 


0.119 


0.106 


Ck 


0.203 


0.122 


0.0371 


0.207 


0.127 


0.0422 


0.211 


0.132 


0.0471 


fiK{GeV) 


0.774 


0.821 


0.869 


0.775 


0.823 


0.870 


0.776 


0.824 


0.871 



TABLE III: Kaon twist-2 wavefunction parameters for some typical Gegenbauer moments, where 
/io = IGeV. Note the obtained WF parameters are for /i = 2.2GeV. 



and 



[l+BKCf\2x-l)+CKC'/'{2x-l)] 



3/2, 



A 



K 



x(l — x) 



exp 



X 



1 — X J 
(34) 

where q = u, d, 2^(1 — 2x) are Gegenbauer poljTiomials. The constitute quark masses 
are set to be: rrig = O.SOGeV and rris = 0.45GeV. After doing the integration over 
the transverse momentum dependence, we obtain the twist-2 kaon DA, e.g. (pxixyfj^b) = 



k_L), where /ib = 2.2 GeV for the present case. The Gegenbauer moments 



twist 


vr 


fx = 1 GeV 


K 


H = 1 GeV 


4 


Si 


0.18 ± 0.06 GeV^ 


Si 


0.20 ± 0.06 GeV^ 






f (0.2 ±0.1) 


ex 


f (0.2 ±0.1) 



TABLE IV: Input parameters for the pion and kaon twist-4 DAs' [1 
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FIG. 2: Typical distribution amplitudes (j)p{x) at = 2.2GeV, where and (pxix) are for 

aJ(lGeV) = 0.115 and aJ(lGeV) = -0.015, of (IGeV) = 0.06 and af (IGeV) = 0.25 respectively, 
and (j)^^"'{x) and (j)^"'{x) are for Gegenbauer expansion ()Aip with the same Gegenbauer moments. 



O'n'^ifJ'b) is defined as 

<'''(/^^) = ^°;t!'T 7.:!.:ii • (35) 



^n,Kr.. ^ _ /o dx(f)^^K{'^ - x, /ib)C^/^(2x 

/o^rfa;6a;(l-a;)[C^/^(2x-l)]2 



The four unknown parameters can be determined by tlie first two Gegenbauer moments, 
the normahzation condition Jq dx Sk\<^il ^T^7r,_ft:(a^5 k^) = 1, and the constraint (k^)^^ ~ 



/ (ixd^kx 


kil 


|*,r,i^(x,kx)P 


/ dxd'^\<i± 


^^,i^(a;,kx)|2 



(ki)y2 = 0.350GeV |29|, where the average value of the transverse momentum square is 
defined as 

(k^ ^^/^ 

Some typical parameters for the pion and kaon WFs are presented in TabHTland Tab lllll A 
comparison with the conventional Gegenbauler expansion DAs is presented in Fig. ([2]). The 
remaining parameters for the twist-4 DA's {^1 ki ^t^,k) are presented in Tab JIVt which are 
taken from [l^ . 

B. Properties of fE'^nxil^) within QCD LCSR with chiral current 

Taking the above mentioned parameters, we discuss the properties of fs'^'^K^'f) within 
QCD LCSRs with chiral current. At the maximum recoil region, = 0, by varying the 
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q (GeV^) 



q(GeVO 



q(GeVO 



FIG. 3: Uncertainties of the S — > vr form factors f^'^'^{q'^) within the allowable regions for the 
undetermined parameters. The center dashed line is for m;, = 4.80 GeV, 02(1^6^) = 0.115, 
al{lGeV) = -0.015, 5l = 0.18 GeV^ and = 0.525. 

parameters within their reasonable regions, we obtain 

/+:^^(0) = 0.260t°:°^^ , fl^M = 0.276lHi (36) 

and 

f^^AO) = 0.334t',Z , f^^AO) = 0.3791°:°??. (37) 

By comparing B —>■ K form factors with the B —>■ n form factors, we find the following 
S'f//(3)-breaking effects among the B — > light form factors: 

/B^ic(O) _ 1 90+0. 06 fs^Ki^) _ -, 07+0.07 foQ\ 

/+,o /n\ ~ -^-^"-o.os ' fT /n^ ~ ^■'^'-om- y'^°) 

It is found that this larger S'f/j(3)-breaking effect is obtained by taking a larger a^(lGeV) G 
[0.10,0.40]; if taking a smaller af{lGeV), then one can obtain a smaller S'[//(3)-breaking 



effect, e.g. = 1.13 ± 0.03 for a^{lGeV) e [0.05,0.10] [sj and = l-OStHj 

for af (IGel/) G [-0.11,0.27] |30|. Note that this, larger S'f//(3)-breaking effect is consistent 
with the some other LCSR calculation as Refs. 
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231 ] and a recently relativistic treatment 



f'*''" (0) 

^^'suii-ouuwmgtji tjquaLiuii in K^y^u, i.e. b~*k 

So a better determination of af' (IGeV) will be helpful to obtain a better understanding of 



that is based on the study of the Dyson-Schwinger equation in QCD, i.e. 43^^ = 1.23 j31 [. 



the S'f//(3)-breaking effect. 

We show the B ^ tt vector, scalar and tensor form factors with their corresponding 
errors in Fig.Q, where the center dashed line is for = 4.80 GeV, a2{lGeV) = 0.115, 
alilGeV) = -0.015, §1 = 0.18 GeV^ and = 0.525. For fB:^^{q^), the lower edge of 
the shaded band is obtained by setting nib = 4.75 GeV, a^{lGeV) = 0.0, al{lGeV) = 0.0, 
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q'(GeV') q^(GeV') q'(GeV') 



FIG. 4: Different parts' contributions to the B ^ tt form factors /^l^'^ {q'^) for all the parameters 
taken to be their center values. The curve of asymptotic total in the middle figure stands for the 
LO /^-+7r(9^) up to twist-3 that is derived from Eg. ([29]) . 

Si = 0.12 GeV^ and e^^ = 0.2625 and the upper edge is obtained by setting mi, = 4.85 
GeV, a^ilGeV) = 0.230, al{lGeV) = -0.030, 6^ = 0.24 GeV^ and = 0.7875. While 
for fB^T,{q'^), the lower edge of the shaded band is obtained by setting rub = 4.75 GeV, 
a^{lGeV) = 0.0, al{lGeV) = 0, 5l = 0.24 GeV^ and = 0.7875 and the upper edge is 
obtained by setting = 4.85 GeV, al^{lGeV) = 0.230, al{lGeV) = -0.030, 5l = 0.12 
GeV^ and = 0.2625. This difference is caused by the fact that the twist-4 structures lead 
to positive and negative contributions to the fB^T^iQ^) and f^^^{q^) respectively. The main 
uncertainties of the form factors are caused by the value of and and it can be found 
that all the B ^ n form factors shall increase with the increment of rrih and . Further more, 
we obtain /J_^^(0)//^_^^(0) G [1.03, 1.08]. This shows that the NLO correction shall affect 
the usual simple relation ( 128|) . e.g. [/B-^7r(0)//^_^^(0)] = (m^ + m^)/m5 G [1.12, 1.14], to a 
certain degree. Furthermore, we show contributions to the B -n form factors fs'^^ni.'f') 
from the different parts in Fig.(j4]), where all the parameters are taken to be their center 
values. For /^^^(g^), it can be found that the LO twist-2, the NLO twist-2 and the LO 
twist-4 contributions are positive, more specifically at = 0, they are about 68%, 26% and 
6% respectively, /^^.^^(g^) is very close to the asymptotic LO result derived from Eq.( l29i) . 
which is caused by the fact that the LO twist-2 gives zero contribution to the sum of the form 
factor [fs^T, + fE^n] a-nd then [/^^^ + /b_,^] gives negligible contribution to /^^^(g^). For 
the LO twist-2, the NLO twist-2 give positive contribution while the LO twist-4 
gives negative contribution, more specifically at = 0, they are about 72%, 30% and —2% 
respectively. 
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q'(GeV') 



FIG. 5: Uncertainties of the B ^ K form factors f^'^Tr io^) within the ahowable regions for 
the undetermined parameters. The center dashed hne is for rrif, = 4.80 GeV, af"(lGey) = 0.06, 
a^{lGeV) = 0.25, 6\ = 0.20 GeV^ and eK = 0.525. 




FIG. 6: Different parts' contributions to the B ^ K form factors fs'^'Kil'^) parameters 
taken to be their center values. The curve of asymptotic total in the middle figure stands for the 
LO f%^j^{q^) up to twist-3 that is derived from Eg. ([29]). 

Second, we show the B K form factors with their corresponding errors in Fig. ([5]), where 
the center dashed line is for = 4.80 GeV, af{lGeV) = 0.06, af (IGel/) = 0.25, Sj^ = 0.20 
GeV^ and ex = 0.525. For f^fij^{q'^), the lower edge of the shaded band is obtained 
by setting = 4.75 GeV, a^{lGeV) = 0.09, a^\lGeV) = 0.10, 5% = 0.14 GeV^ and 
ex = 0.2625 and the upper edge is obtained by setting nib = 4.85 GeV, af{lGeV) = 0.03, 
a^ilGeV) = 0.40, 5| = 0.26 GeV^ and = 0.7875. While for /^^^(g^), the lower edge of 
the shaded band is obtained by setting mfe = 4.75 GeV, af{lGeV) = 0.09, a^{lGeV) = 0.10, 
6]^ = 0.26 GeV^ and ex = 0.7875 and the upper edge is obtained by setting rrih = 4.85 
GeV, af (IGeV) = 0.03, a^{lGeV) = 0.40, 6l = 0.14 GeV^ and ex = 0.2625. The main 
uncertainties are caused by the value of m^, af^ and af^, and it can be found that all the 
B ^ K form factors shall increase with the increment of and af^, and decrease with 
the increment of af . As for the LO results, we obtain [/B-^/^(0)//^_^^(0)]i^o ^ [1-19, 1.22]; 
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while for the NLO results, we obtain [/J_^/^(0)//^_^^(0)]ArLo ^ [1.09,1.15]. Furthermore, 
we show the different parts' contributions to the B K form factors /^^'^(g^) in Fig. ([6]), 
where all the parameters are taken to be their center values. For fs'^K^'f')^ can be found 
that the LO twist-2, the NLO twist-2 and the LO twist-4 contributions are positive, more 
specifically at = 0, they are about 67%, 27% and 6% respectively. Even though the 
LO twist-2 gives zero contribution to the sum of the form factor [fs^x + /b-^a'] but due 
to S'f//(3)-breaking effect they shall give sizable contribution to fB^xiq"^), so fs^Kio^) 
higher than the LO result derived from Eg. (129!) as shown in Fig. ([6]). For f'^_^f^{q'^), the 
LO twist-2, the NLO twist-2 give positive contribution while the LO twist-4 gives negative 
contribution, more specifically at = 0, they are about 70%, 32% and —2% respectively. 

IV. COMPARATIVE STUDIES OF /^I^J^lg^) WITH OTHER APPROACHES IN 
QCD LCSRS 

A. A striking advantage of the present approach with the chiral current 

The adopted chiral current approach has a striking advantage that the twist-3 LC func- 
tions which are not known as well as the twist-2 light-cone functions are eliminated, and 
then it is considered to provide results with less uncertainties. On the other hand, by using 
the standard weak current in the correlator as shown by Eqs. lfTpj) . it has been pointed out 



that the twist-3 contributions can contribute ~ 30 — 40% to the total contribution [32] . 
So to obtain a more accurate result, one has to calculate the above correlator by including 
one-loop radiative corrections to both the twist-2 and the twist-3 contributions. Such a 
calculation together with the updated pion and kaon twist-3 wave functions has been done 
by Ref.fl. 

It may be interesting to do a comparison of their results with our present ones so as to 
show whether these two treatments are consistent with each other or not. For such purpose, 
we adopt the following convenient form for the QCD sum rules obtained by Ref. 5{], which 
splits the B ^ P form factors into contributions from different Gegenbauer moments: 



B^P 



= r{q') + af (/^o)/<(g') + «r(/io)r"(g') + ar(/xo)/<(g'), (39) 



where contains the contributions to the form factor from the asymptotic DA and all 

P P P 

higher-twist effects from three-particle quark-quark-gluon matrix elements, /"^i '"^a '"^^ con- 
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FIG. 7: f^'^'^iq^) with the allowed values of and a J being correlated and given by the rhomboid 
shown in Fig{TJ The solid line is obtained with 02 (IGeF) = 0.0 and a\[\GeV) = 0.0 and the dashed 
line is obtained with 02 (IGey) = 0.23 and a\{lGeV) = —0.030, which set the upper and the lower 
ranges of fs^'nil'^) I'espectively. As a comparison, the shaded band shows the results of Ref.^] 
together with its 12% theoretical uncertainty. 




q'(GeV') q'CGeV") q'(GeV') 



FIG. 8: fs'^Kiq^) for a^{lGeV) G [0.03,0.09] and af (IGeF) G [0.10,0.40]. The solid line is 
obtained with af{lGeV) = 0.09 and a2{lGeV) = 0.10 and the dashed line is obtained with 
a^{lGeV) = 0.03 and a2 (IGeV) = 0.40, which set the upper and the lower ranges of /^^'^(g^) 
respectively. As a comparison, the shaded band shows the results of Ref.jsl together with its 15% 
theoretical uncertainty. 

tains the contribution from the higher Gegenbauer term of DA that is proportional to af , 
af and af respectively. The explicit expressions of J'^^-'^f -"f-'^f for all the mentioned form 
factors can be found in Table V and Table IX of Ref. [sl. And in doing the comparison, we 
take the same DA moments for both methods. 

We show a comparison of our results of with those of Eq. fl39|) in Figs. fl71l8|) 

respectively. Fig.(l7j) shows /^^'J(g^) with and being correlated and given by the 
rhomboid shown in Fig.([l]), where the solid line is obtained with a2{lGeV) = 0.0 and 
al{lGeV) = 0.0 and the dashed line is obtained with a^{lGeV) = 0.23 and al{lGeV) = 
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—0.030, which set the upper and the lower ranges of respectively. Fig.® shows 

ft'^l{q^) with a^{lGeV) G [0.03,0.09] and a^{lGeV) G [0.10,0.40], where the solid line 
is obtained with a^ilGeV) = 0.03 and af^(lG'e\^) = 0.10 and the dashed line is obtained 
with af^{lGeV) = 0.09 and af"(lGeV^) = 0.40, which set the upper and the lower ranges of 
fs'^Ki.'f') respectively. As a comparison, the shaded bands in these figures show the results 
of Eq.dnnD within the same <x^ and region and with their estimated [12% + 3%] theoretical 
uncertainty, where the extra 3% uncertainty is from uncertainty |5|. 
More explicitly, we show the comparison in detail: 

• At the large recoil region = 0, Ref.Q gives fst^i^) = 0-258 ± 0.031, /^_^^(0) = 
0.253 ±0.028, /b£a'(0) = 0.304 ±0.076 and fB^AO) = 0.332 ±0.080. It can be found 
that our results as shown by Eqs. fl36f37l) are consistent with those of Ref. [5|], especially 
in the lower region. 

• With the increment of g^, the form factors of Ref. 5| increase faster than ours. We can 
see this clearly from the scalar and tensor form factors fs'^nKil'^)- These differences, 
especially in the larger region are mainly caused by the treatment of the twist-3 
contribution and by the different treatment of the uncertainty. The twist-3 contribu- 
tion can affect the shape of the form factors. For example, as shown by Fig.(j4]), the 
present obtained /^^^^(g^) are close to the LO result derived from Eq. fl29l) : while 
Ref. jsi gives a larger f%^^j^{q^) at higher g^ region due to the fact that the twist-3 
contribution to [fB^T,j<: + Ib^tt^k] is dominant over the leading twist contribution at 
large momentum transfer. In Ref.|^ the total uncertainty is obtained by adding up 
the uncertainties caused by each parameter in quadrature; while at the present, we 
vary the parameters within their possible regions and adopt the minimum and the 
maximum ones as the uncertainty boundary. Moreover, we have adopted a simple 
overall uncertainty 12% or 15% for the form factors within all g^ for the LCSRs of 
Ref. jsj, which in fact should be varied according to different g^, e.g. we have found 
that such uncertainty may be up to 5% for the mentioned form factors with the region 
of g2 e [0, 14] GeV^. 

• One may observe that in the lower g^ region, different from Ref.j^ where F^~Q'^[q^) 
increases with the increment of both af and af^, our present predicted F^~^:^{q^) will 
increase with the increment of af' but with the decrement of . This difference is 
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caused by the fact that we adopt the pion and kaon DAs derived from their wave 
functions to do our discussion, whose parameters are determined by the combined 
effects of and af-; while in Ref.jsl, af" and af^ are varied independently and then 
their contributions are changed separately. 



B. A comparison of the choosing of pole or MS 6-quark mass 



Refs. Q, [l7| has argued to used MS 6-quark mass instead of the pole quark mass. The 
MS 6-quark running mass (mf,) is related to the one-loop 6-quark pole mass {ml) through 
the following well-known relation: 



n 



*2' 
b 



m,(/i) = <|l + ^^^|^-4 + 31n-^j|. (40) 

With the help of the relation pUj) . one can conveniently transform the form factor expressions 
among these two choices of 6-quark mass. And one only need to be careful to use all the 
parameters calculated under the same choice, e.g. the value of should be calculated by 
using the same currents in the correlator and under the same choice of 6-quark mass. By 



calculating the ordinary correlators ffT|2]) up to NLO and by varying the MS 6-quark mass 



within the region of mfe(mfe) = 4.164 ± 0.025 GeV, Refs.pJ, IjJ obtain fstA'^) = 0-26t[j 



h0.04 
1.035 



fl^M = 0.255 ± 0.035, fttA^) = 0.36t°:[Ji /|^^(0) = 0.38 ± 0.05, and 

fp^K^^) -I Qo+0.11 fp^Ki^) 1 /IO+0-18 (■A-i\ 

These results are consistent with ours and also with those of Ref. within reasonable errors, 
which is also calculated by taking the pole quark mass. This shows that these two choices 
of 6-quark mass are equivalent to each other. 



C. Extrapolations of the LCSR results to higher region 

In order to allow a simple implementation of our results, we present a parametrization 
that includes the main features of the analytical properties of the form factors and is valid 
in the full physical regime < < (m^ — mp)^. Following the same argument of Ref.j^, 
we fit the LCSR results to the following parametrizations that are based on the procedure 



advocated by Becirevic and Kaidalov [33|, where we take the LCSR results with all the 
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parameters taken to be their center values to do the extrapolation, i.e. the 6-quark one-loop 
pole mass = 4.8 GeV, a^{lGeV) = 0.115, al{lGeV) = -0.015, af{lGeV) = 0.06 and 
af^(lG'eV^) = 0.25. To measure the quality of the fit, we introduce the parameter A that is 
defined as 

fit) - f\t) 



A 



100 max 

t 



fit) 



1 2S 
te{0,-,...,-,12\ GeV^ 



(42) 



i.e. it gives, in per cent, the maximum deviation of the fitted formfactors from the original 
LCSR result for < 12 GeVl 



for fly. 



fi<f) 



ri 



r2 



I {m\) 



1 



(43) 



fit 



where = 5.325 GeV 18|] is the mass of i?*(l ). For the fit parameters are 
ri = 0.7411, ra = -0.4815 and ml^ = 40.01 GeV^ for A ~ 0.05. And for the fit 
parameters are ri = 0.7742, r2 = -0.4952 and ml^ = 34.71 GeV^ for A ~ 0.9. 



for /^^: 



fiq' 



l-gV(mf)2 



+ 



r2 



1 - gVm|t' 



(44) 



181 ] is the mass of the -8^(1 ). For f^, the fit parameters are 



where mf = 5.413 GeV 
ri = 0.8182, r2 = -0.4862 and m|, = 41.61 GeV^ for A ~ 1.3. And for /f , the fit 
parameters are ri = 0.893, r2 = —0.5073 and m^^ = 33.13 GeV^ for A ~ 1.8. 



for /o 



foiq' 



r2 



1 — /m? 



fit 



(45) 

46.09 GeV^ for A = 0.07. For 



For the case of pion, we obtain r2 = 0.2596 and m^^ 
the case of kaon, we obtain r2 = 0.332 and m|t = 61.64 GeV^ for A = 1.3. 



A comparison of the lattice QCD results can 
results in the literature for B 



vr, e.g. 



34 



35 



3e found in Fi g.(p ). There are many Lattice 
40| and etc., for convenience. 
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we have taken the unquenched lattice QCD result f36] and the quenched lattice QCD result 
351] oi B ^ -K form factors. While for the B —>■ K form factors there is little lattice QCD 
results, and we adopt the prehminary results derived by Ref.[4(| to do our discussion. 
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FIG. 9: Extrapolations of the LCSR Results of i? ^ vr and K form factors for higher . For com 



parison, Left diagram shows the unquenched lattice QCD result 
lattice QCD result 



3q | (Diamond) and the quenched 



35l | (Triangle) with their corresponding errors for the vector and scalar i? — > tt 

4^ 



form factors; Right diagram shows the Lattice QCD results 
form factors (Asterisk). 

V. SUMMARY 



for the vector and scalar B ^ K 



We have calculated all the B ^ P transition form factors F^f.piq'^) with chiral current 
in the LCSR up to NLO, in which the most uncertain twist-3 contributions have been elimi- 
nated naturally, with the b-quark pole mass that is universal. The S'f//(3)-breaking effects in 
B ^ K form factors have been carefully discussed and their values depend on the moment 
of the kaon distribution amplitude, af^(lG'ey). It is found that = l-28to;o8 ^^"i 

= 1.37iH2 for af{lGeV) E [0.03,0.09] and a^{lGeV) E [0.10,0.40]. Based on the 
LCSR with chiral current, we have made a comparative study on the properties of transition 



form factors with those obtained in literature p, 14, 17], in which the radiative corrections 
on both the twist-2 and twist-3 parts should be treated in equal footing. It has been found 
that the present results are less uncertain under the same parameter regions to consider the 
radiative corrections since the twist-3 contributions have been eliminated naturally in the 
adopted method, so our results are simpler and consistent with those in literature that have 
been derived with the usual correlators. These form factors are important ingredients in the 
analysis of semileptonic B decays, especially our present results may be helpful to clarify 
the present conditions for the B 'r]^'\i^i'i, i^i") decays and then a better understanding 
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of the 7] and t]' mixing 



m 
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APPENDIX A: PION AND KAON DISTRIBUTION AMPLITUDES 

Generally, the pion and kaon twist-2 and twist-4 DAs can be written in the following 
forms: 

• twist-2 DAs: 



ifp{u,^) = 6uu l + a'i'{fi)C^^\2u-l) + a^{fi)C2^\2u-l)+a^{^)Ct'\2u~l) + 



3/2/ 



(Al) 



where P stands for tt or K respectively, ■ ■ ■ stands for even higher Gegenbauer terms. 
twist-4 DA's Q: 



^4p(ai) = 30a3(Q;2 — cti] 



-30a^ 



hoo + hoia-s + ^ /iio(5a3 - 3) 
^oo(l - «3) + V "3(1 - tta) - Qaia2 



+hio^a3{l - as) - ^{aj + a^^ 

$4p(aj) = 120aia2a3 [aio("i - "2)] , 
^Ap{oii) = 120aia2a3 boo + fio(3a3 - 1)] 



(A2) 



where 



^00 

^^10 

hoi 



voo = - 
21Ml 



Ml 



51 



-rji^pLOAP 



21M| 



9 

20 

s:2 



Tji^pUJiip — Sp€p 

3 



-r]/ipLOiP 



20 



20 



20 



23 



and 



no 



-^ViP^4P + Mp = -dpep + —a^ Mp, 



with = 6p/Mp, U4P = 8ep/21. Taking the leading meson-mass effect into consid- 
eration, the remaining two-particle DA's of twist 4 can be written as Sj: 

. , 4:UU 



Hp[u 



u 



g I — 5uu [30/ioo + 4/?.oi(3 + uu) + 5/iio(— 3 + 2mm)] 
+2aio[6 + mm(9 + 40mm)]} + 8aio{2M^(10 - 15m + Gm^) In 
+2m3(10- 15M + 6M2)lnM}, (A3) 
5 - 4/ioo - 2/ioi + /iio + 4 (-4aio + 6/ioo + 4/ioi + ho) u 
-6 (4/ioo + Qhoi + 9/iio - 16aio) + 20 (2/ioi + 5/iio - 8aio) m^ 
+10 (Saio - 2 V - 5/iio) M^l . (A4) 



Setting Mp and Mp ^ M/^, one can obtain the pionic and the kaonic twist-4 
DAs respectively. 
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